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Abstract. Jets and outflows are known to form out of advective accretion flows around 
black holes. Hard photons from the centrifugal barrier directly hit the electrons and 
deposit momentum on them. For optically thick flows such deposition is not efficient, 
but for optically thin flows matter could be accelerated to relativistic speed. In fact, 
even bound matter could be made free through successive deposition. We discuss these 
possibilities. 



Proceedings of Heidelberg International Gamma- Ray Astronomy Conference (Eds. 
F. Aharonian and H. Voelk) 

INTRODUCTION 

The effect of radiative momentum deposition on astrophysical outflows around 
compact objects were studied using early disc models by many workers (e.g. Vincent 
Icke [1] and Abramowicz and Piran [2]). Chattopadhyay and Chakrabarti [3,4] , 
revisited the problem using the radiation from advective discs [5] . The outflow was 
found to be accelerated appreciably, pushing bound matters to infinity. Compact 
objects do not have atmospheres like stars do, therefore the outflows or jets have 
to be generated from the inflowing matter itself. As accreting matter comes closer 
to the compact object, it is slowed down by the centrifugal force and may cause 
shock [6], forcing the flow to jump discontinuously from supersonic to subsonic 
state. This hot, slowed down region puffs up in the form of a torus, known as the 
centrifugal pressure supported boundary layer (hereafter, CENBOL). Chakrabarti 
and Titarchuk [5] pointed out that radiation from this region may dictate whether 
a black hole is in a soft state or in a hard state. Similarly, Chakrabarti and co- 
workers [7,8] pointed out that the same region also causes formation of outflow. 
In this region, the subsonic outflowing matter is continuously bombarded by hot 
photons and hence apart from usual thermal acceleration, radiative acceleration 
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Fig. 1: Schematic diagram of CENBOL and disk/outflow geometry. D is 
the source point, B is the field point. O is the position of compact object. 
The local unit normal h is along DC. OO =Z ,0 D = X 



should be important too. Thus the toroidal region or the CENBOL of accretion 
disc not only emits high energy radiation, but its typical shape helps the radiation to 
be focussed on the axis of symmetry. It is of no surprise that we find an impressive 
degree of acceleration at this region. Our present exact computation by radiative 
momentum deposition (RAMOD) also vindicates earlier investigation [3,4] , that 
momentum deposition force can also drive bound matters as outflows. 



Solution of two-temperature equations inside a centrifugal barrier supported shock 
located at around r = 10r s (lr p = 2G jf B ) is presented in Chakrabarti and 
Titarchuk [5]. The inner surface of the CENBOL can be approximated by a cone, 
which we do and compute RAMOD along the axis of symmetry. 
Figure 1, gives a schematic diagram of the disc/outflow geometry. We choose cylin- 
drical (r,0,z) coordinates, where the compact object is at the origin O. We assume 
the entire disc/outflow system is symmetrical about the axis of rotation of the disc, 
therefore RAMOD has only two components (D r and D z ). We define a <fi = plane 
PMNO The flux of radiation at field point B(r,z)i on PMNO, from differential 
area dA around source point D(X, 0, Z') is given by: dF = 1 cos L(EDB)dVt ) where 
/ is the frequency averaged radiation intensity on the surface and dQ is the solid 
angle subtended by the differential area dA at B and is equal to dA/BD 2 . The 
RAMOD will be along the flux i.e, along BD and will have two components, one 
along the projection of BD i.e, AB and the other transverse to it. The transverse 
component will cancel each other as we integrate on <fi. Hence the two components 
of RAMOD is given by; 



CALCULATION OF RAMOD 



D z (r,z) = D J j[ 
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and, 



D r (r,z) = D J J[ 



BD 2 + ED 2 - BE 2 , X „ AG 



d(f)dX 
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where, D = ^ £2: , 9 =semi vertical angle of the inner surface of CENBOL, 

' u cm p ' & ' 

m p =mass of the proton, c = velocity of light. To estimate / we assume that 
the dissipation of binding energy of the last stable orbit of accreting matter is the 
source of radiation. Hence, / = q^-j, where L = i]M acc c 2 , is the amount of rest 
mass energy of the accreting matter converted into radiation per unit time. Here 
M aC c is the accretion rate, rj is the convertion ratio and is equal to 0.06 for the 
Schwarzschild metric. Qt = 27r, is the solid angle in which the radiation is coming 
out locally. A = 2tiR(H 2 + R 2 ) 1 / 2 , is the total of conical CENBOL (see, 
Fig. 1). In our calculation of RAMOD we do not consider the attenuation of the 
radiation field. Optically thick flow would be difficult to accelerate this way. Our 
assumption of optically thin flow (i.e., r < 1) puts limits on the mass outflow rate 
and the mass of the central object. The density of the radial outflow is given by; 

M , . 

P = — , (3a) 

where r =radius vector of outflow in spherical geometry, v =radial velocity. Let 
r = L = 10r g be the length scale where the above approximation is valid. We 
rewrite the expression of the density using the definition of sound speed, 

In a radiation dominated flow, P gas « P raa - i.e, 

aT 3 jj,m p 

P<< ^ JL ' (4) 

which provides a upper limit of density. Combining expressions (3b) and (4), we 
get; 



M 6 an 4 m p 4 
L 12 ^k B 4 
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(5) 



The condition for the optically thin flow (r < 1) in the Thompson scattering regime 
for length scale ~10r s , also provides the estimate of upper limit of density, 

1.69xl0 27 o 
P < ^ gm cm ' . (6) 

From (5) and (6), in order to have optically thin outflows, we get: 

M 6 



< 1.7xl0" 70 gm sec~ d (7) 
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or, 



m w m w < 6.1 xlO" 8 , 



(8) 



Fig. 2: The vector field of radiative force due to the radiation from CEN- 
BOL. M aC c = IOMedd- The position of shock 10r g , specific angular mo- 
mentum of the accreting matter 1.6 ^—^ , H = 4.4r 9 , R = 6.671r g . 

where, rhio is the mass outflow rate (M) in the units of Eddington rate for a 
1OM black hole and mio is the black hole mass (Mb) in units of 1OM . Increas- 
ing M acc increases disk luminosity and hence magnitude of RAMOD. However, 
Chakrabarti [9] showed that higher M acc results in higher outflow rate M, which 
results in outflows with higher density p and therefore higher r. In case of optically 
thick outflows radiative acceleration is not efficient, as the intensity of radiation 
from the disc will be attenuated by a factor e~ T while passing through the flow. 
In our work the outflow considered are optically thin and hence respects equa- 
tions (7) and (8). On integrating equations (1) and (2) from <p = 0—>-2tt and 
X = l.5r g —*R, we get the value of RAMOD everywhere. Figure 2, shows the 
vector field of RAMOD. We choose M acc = IOMedd- We have also chosen the 
shock location for accretion to be at a distance 10r 9 [5] on the equatorial plane of 

the accretion disc, the specific angular momentum of the accreting matter to be 
I qZGMb^ and thug we have H = AAr ^ R = 6 6 7 lr . g 

From Fig. 2 we see that very close to the horizon the radiation force is directed 
towards it. As one goes away, RAMOD becomes positive it is still directed towards 
the axis of symmetry. Further out it falls off ~l/r 2 and force field becomes more 
and more radial. This focusing effect, will tend to collimate the outflow. 

GLOBAL SOLUTION WITH RADIATIVE MOMENTUM 

DEPOSITION. 

As a special case of solutions of equations (1) and (2), RAMOD along the axis 
of symmetry i.e. D z (0,z), was calculated by Chattopadhyay and Chakrabarti [4] 
and was approximated as D r ~D z (0, z), for conical outflow very close to the axis of 
symmetry. 



V r = J D [cos2($ 1 - 6) - cos2($ 2 - 6)} 



(9) 
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Fig.3(a-b):Variation with log(r) of (a) v, RAMOD corresponding to 
WMedd (solid) and Bondi outflow (long dashed); (b)Mach number M, 
RAMOD corresponds to M acc = 15Medd (solid) and thermally driven 
outflow (long-dashed). RAMOD solution becomes unbound. 



where, $ x = tan 1 [^^^] and, $ 2 = tan 1 [ 2(1 r ^ 0) 
H is the maximum height of the CENBOL from the equatorial plane and, D = 
f^ 21 , where J, a-r, c, m v has there usual meaning. Using expression (9) in the 
momentum balance equation, and in the entropy generation equation we solve for 
flow topology of the radial outflow. Figure 3a shows the comparison of radial veloc- 
ity v with and without the RAMOD. Velocity is plotted against log(r). The solid 
line corresponds to IOMedd and the long dashed line is the solution of Bondi [10] 
outflow, i.e. outflow not acted on by RAMOD. The initial velocity at which flow is 
launched for the two cases is 0.1717c at radial distance 2.598r 5 . The v t for the above 
cases are 0.743c, 0.5804c respectively. We see that RAMOD accelerates the outflow 
and higher RAMOD results higher terminal speed. RAMOD not only accelerates 
but also increases the energy of the flow often freeing bound matter. Figure 3b 
shows the comparison of Mach number M{= v/a) variation variation with log(r). 
The long-dashed curve represents the thermally driven negative energy flow, and 
the solid curve represent the outflow which is acted on by RAMOD correspond- 
ing to M acc = 15Medd- The initial parameter for both is, Ei n = —0.0349c 2 at 
ri n = 2.509r 9 . In this case of later, RAMOD increases the energy of the flow to 
positive value, pushing it to infinity as free outflow and making it supersonic in the 
process, while the former remains bound and ultimately dives back to the compact 
object. 

This work is partly supported by DST project Analytical and Numerical Studies 
of Astrophysical Flows Around Compact Objects. 
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